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Abstract 

We present calculations of magnetic exchange interactions and critical temperature Tc in 
Gai_a;Mn^As, Gai-xCr^^As and Gai-^Cr^^N. The local spin density approximation is combined 
with a linear-response technique to map the magnetic energy onto a Heisenberg hamiltonion, but 
no significant further approximations are made. Special quasi-random structures in large unit cells 
are used to accurately model the disorder. Tc is computed using both a spin-dynamics approach 
and the cluster variation method developed for the classical Heisenberg model. 

We show the following: (i) configurational disorder results in large dispersions in the pairwise 
exchange interactions; (ii) the disorder strongly reduces T^, (Hi) clustering in the magnetic atoms, 
whose tendency is predicted from total-energy considerations, further reduces Tc. Additionally the 
exchange interactions J{R) are found to decay exponentially with distance on average; and 
the mean- field approximation is found to be a very poor predictor of Tc, particularly when J{R) 
decays rapidly. Finally the effect of spin-orbit coupling on Tc is considered. With all these factors 
taken into account, Tc is reasonably predicted by the local spin-density approximation in MnGaAs 
without the need to invoke compensation by donor impurities. 
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Dilute magnetic semiconductors (DMS), i.e. semiconductors doped with low concentra- 
tions of magnetic impurities (usually Cr, Mri, or Co), have attracted much interest because 
of their potential application to spintronicsjl|, 0] . Gai-^Mn^-As is the most widely stud- 
ied DMS, and it continues to attract interest because it is one of the few DMS where it is 
generally agreed that the magnetism is carrier-mediated. (This is important in spintronics 
because the magnetic state can be manipulated by electrical or optical means.) 

In recent years Curie temperatures in Gai-ajMn^jAs have risen steadilj^ reaching ~170K 
for x~0.08 when grown in thin films annealed at low temperature 0, S]. It is generally 
believed defects (probably Mn interstitials) migrate out of the as-deposited films during the 
anneal, largely eliminating donor defects that hamper ferromagnetism. Since most practical 
applications of spintronics require room-temperature operation, a crucial question is then, 
what is the ultimate limit to Tc in the DMS compounds, and in Gai-xMn^^As in particular? 

This question was first addressed by Dietl in his now classic paper [oj, where he predicted 
a wide range of Tc in tetrahedrally coordinated alloys. This stimulated a great deal of 
interest, although there is a growing consensus that most of the claims of thatpaper were 
artifacts of the assumptions in his original model. On the other hand, Akai 3] first used 
the local spin-density approximation (LSDA) to estimate within the Coherent Potential 
Approximation (CPA) in (In,Mn)As; he argued that a double exchange mechanism was a 
more appropriate description of the magnetism than the pd exchange assumed by Dietl. 
Since then 
of groups 

energies of a fixed atomic but multiple-spin configurations, or by a linear-response technique 
within the CPA. 

To date, disorder has almost always been neglected or treated within some mean-field 
(MF) approximation (MFA), either in the computation of the exchange parameters them- 
selves, or in the subsequent analysis of magnetization M(T) at finite-temperature, or both 
(though better treatments within k ■ p theory has been reported ^^). The LSDA-I-MF 
predict a rather high Tc for Gai^^-Mn^As (typically 350~400 K for x~0.08 [3|)- The large 
discrepancy with experiment (at least in Mn:GaAs) is usually attributed to the very large 
numbers of compensating defects in real samples, which reduce Tc 3,0,0]. The situation re- 
mains somewhat uncertain because the number of defects still remaining in the best samples 
to date is not known. 



jSDA calculations of exchange interactions have been performed by a variety 

mum 

, usually extracting exchange parameters by calculating total 
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This Letter addresses the issue of the ultimate hmit to Tc in some DMS alloys (focusing 
on Mn:GaAs) by adopting relatively rigorous approach to the calculation of the magnetic 
exchange interactions and T^. Random alloys are approximated by large (128-250 atom) 
supercells where special quasirandom structures (SQS) jl^ are used for the cation sublat- 
tice. Using a 



method 
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inear-response technique within the LSDA and the linear-muffin-tin orbitals 
17|, the magnetic energy is mapped onto a Heisenberg formjisj] 

H = -J2JiR,,)e,-e, (1) 
where the sum is over all pairs ij of magnetic atoms. To model M{T) and T^, Eq.([Q) is 



treated classically and integrated using a spin-dynamics ( SD) techniq ue |l9l| : alternatively 
M{T) is estimated by the cluster variation method (CVM)j20| adapted j2l| to solve Eq.(IT]). 
Thus it is evaluated without recourse to empirical parameters or to the MFA. We show that 
the widely used MFA turns out be a very poor predictor of M(T) in these disorded, dilute 
alloys, dramatically overestimating Tc. 

With SQS we can rather precisely mimic a fully random configuration, but it is also 
possible to consider configurations that deviate from random. This can be important because 
LSDA predicts a strong attractive interaction between magnetic elements 
a tendency towards clustering. In brief, we show that 



which implies 



• the disorder induces large fiuctuations in = J{Rij) for every connecting vector Rij] 

• The fiuctuations in Jij reduce Tc relative to the configurationally averaged Jij = {Jij); 

• clustering reduces Tc, while ordering of the magnetic elements increases Tc. 

Fig. C] shows Jij computed for an ensemble of 108-cation (216-atom) random supercells 



following the method of Ref. 17[, for Gai-^Mn^-As and Gai^^-Cr^N alloys at x=A.Q% and 
x=8.3%. 3x3x3 /c-points were used, enabling the calculation of J to very distant neighbors. 
We chose these two alloys because they are approximately representative of limiting cases. 
For Cr:GaN, the GaN host has a wide bandgap, and the Cr t2 level falls near midgap. It 
broadens into an impurity band with 1/3 occupancy, and is believed to be responsible for 
the ferromagnetic exchange. For Mn:GaAs, most of the weight of the Mn t2-derived state 
falls below the valence band maximum. A second t2 impurity band about 0.1 eV above 
the valence band maximum is mainly responsible for the ferromagnetic exchange coupling 



0.1 

0.01 
0.01 

0.1 



— 1 1 1 — 

Mn:GaAs • x=.083 
* x=.046 



i • • I i ?>ii 

• I • • ■ r 

: I I % • ' 




1.0 1.5 2.0 
Rij/a 



2.5 




FIG. 1: Pair exchange interactions J{Rij), in mRy, for Mn:GaAs and Cr:GaN at two different 
concentrations as a function of Rf,. R,:^ is measured in units of the lattice constant a. 



in this case; the strength of J{R) depends critically on the amount of Mn character in this 
bandQ]. Katayama-Yoshida used the x-dependence of Jq = J2r J{R) (computed within the 
CPA) to identify the ferromagnetism obtained from LSDA with model theories [l^^. Within 
the CPA, Jo ~ 3;^/^ for Cr:GaN, which corresponds to a double-exchange model, while 
Mn:GaAs displays character intermediate between Jq ~ x^^^ and the pd exchange (Jq ~ x) 



usually assumed hj k ■ p models j6 
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Comparing Cr:GaN to Mn:GaAs, Cr:GaN shows substantially stronger nearest-neighbor 
(NN) interactions, owing to its small lattice constant; however J{Rij) decays much more 
rapidly with Rij. This is because the wave function overlap between transition metal d 
states decays much more rapidly for midgap states than near band-edge states. Evident 
also is the large dispersion in Jij for fixed Rij (note J is drawn on a log scale): the root- 
mean square fluctuations AJij = J (jf, — M) are roughly comparable to J. However A J, 



increases with x, and is substantially larger for the wide-gap case (Cr:GaN). Note that there 
is little evidence in either Cr:GaN or Mn:GaAs for oscillatory RKKY-like behavior, which in 
the simplest approximation predicts J{R) ~ cos{2kFR)/R^. Instead, J{R) decays roughly 
exponentially in R^, corresponding to a Fermi surface with imaginary wave number, as would 



24|. 



obtain if the coupling were described by tunneling via a disordered impurity band 

We now apply Eq.(^ to compute M(T), focusing on Tg. Mean-field theory, which es- 
timates the effective field at each site from the average field contributed by other sites, 
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predicts Tc well above room temperature both in Mn:GaAs and Cr:GaN|l3|, |25|. In spite of 



the rather strong differences in the form of J{R) (Fig. ^ 
Mn:GaAs and Cr:GaN have roughly similar for x~0.08 



mean-field theory predicts that 
2511 . This is because the NN inter- 



action in the 



after case is strongest, but the J decays faster with R, leading to a comparable 



mean-field j26| estimate T^^^^ . 

But it should be evident from Fig. ^that the MFA is of questionable reliability. First, 
it is well known that for dilute alloys there is a percolation threshold for the onset of 
ferromagnetism. (The threshold in the present case cannot be readily mapped to known 
models because J{R) is nonneglible for a rather large number of neighbors.) Moreover, the 
large fluctuations AJ{R) may strongly affect T^, especially since AJ{R) itself is purely a 



function of t.e env.onn,en.a, and consequent, of tKe >oca> pe.co.at.on pat. 

To obtain a precise estimate for M(T) and T^, we adopt a spin-dynamics approach [igj]. 
A 200 atom SQS structure (250 atom for the 4% alloy) was used to mimic the random alloy. 
From the TM atoms in the SQS structure, a supercell containing ~ 2000 Mn or Cr atoms 
was constructed to make a simulation cell for prosecuting spin-dynamical simulations. Fol- 
lowing the method described in Ref.^^, the Landau-Lifshitz (L-L) equation was integrated 
numerically at a fixed temperature allowing the system to equilibrate, followed by a simula- 
tion for ~ 2 X 10^ atomic units. The L-L equations were integrated with the Bulirsch-Stoer 
method. As the L-L equation is a first-order equation, global deamons were used for the heat 
bath^^, to ensure ergodic behavior. The average magnetization M(T) was computed as a 
function of temperature, and Tc was estimated from the inflection point in M{T). Owing to 
finite-size effects and the stochastic character of the simulation, Tc could be determined to 
a precision of ~5%. 

Also we employed a CVM approach recently adapted to the classical Heisenberg 
hamiltonion[2l|]. This relatively simple scheme has been found to be accurate in simple 3d 
magnets, overestimating Tc by ~5% (similar to the usual CVM for the Ising hamiltonion 27[). 
We can check the validity both methods in the DMS case by comparing their predictions 
of Tc. Fig 121 shows Tc determined by both methods for Gai_a.Mn^As and Gai_^.Cr^As : 
agreement between the two methods is ~10%, which is quite satisfactory considering the 
complexity of the Jij. T^J^^ is also shown: evidently the MFA rather badly overestimates 
Tc- T^^^^ > Tc by ~200K in the Mn:GaAs alloy, and by a somewhat larger amount in 
Cr:GaAs. The discrepancy is still more dramatic in Cr:GaN (not shown); we find Tc < 50K 
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for all concentrations studied while T^^^ ~ 600K'25]. Indeed we have found this generally 
to be the case when J{R) decays rapidly or when AJ{R) / J{R) is not small. 

These results stand in stark contrast to the ~15% discrepancy between T^^^ and 
typically found in simple metals. The reason is easily understood by considering the effective 
field a mean-field atom sees, H^^^ = J2j Jtj^^j- From the exponential decay of J{R), it is 
evident that Hi will be dominated by the nearest neighbors. But for dilute alloys, near- 
neighbors are not sufficent to form a percolation path. This is immediately evident in the 
extreme case of a NN pair of magnetic atoms well separated from any other magnetic atoms: 
the contribution to T^^^^ from this pair would be high, even though the pair would actually 
contribute nothing to ferromagnetism. 

In Ref. [l^ a small discrepancy between T^^^ and a more sophisticated calculation for 
Tf. was reported. In that calculation the CPA was used to construct an average Jij and 
M{T) modeled by constructing a fee lattice of magnetic atoms, using concentrated-weighted 
Jij for the exchange parameters. It would seem that their conclusions are an artifact of the 
neglect of configurational disorder (except in the computation of Jij). Better would be to 
estimate Jij within the CPA, and then construct a disordered simulation cell using the Jij 
to estimate M{T). Still this approach neglects fluctuations A J, which as we have seen are 
comparable to Jij itself. To assess the effect of fluctuations, we repeated the calculation 
for Tc within the CVM, replacing the environment-specific Jij with the configurationally 
averaged Jij. For Gai_2;Mna;As at a;=0.08, the effect of disorder {Jij Jij) was to reduce 
Tc by 50 K. (It is interesting that the MFA predicts the opposite trend, because of an artificial 
tendency for M^^^^{T) to track whichever site i has the largest Hi . Then T^^^^ — 



is positive 



26| and increases with A J/ J. This explains why a tight-binding+MF analysis 



predicted that disorder increases Tc.) 

We next consider the effects of nonrandomness. As noted above, real DMS alloys should 
exhibit some clustering owing to the attractive interaction between magnetic elements j^. 
The true situation is complicated by the nonequilibrium growth required to stabilize the alloy 
in the zincblende structure. Nevertheless the Mn-Mn or Cr-Cr binding energy is calculatedjq] 
to be an order of magnitude larger than the growth temperature (~250K), and some pairing 
or other clustering should be expected, particularly since films must be annealed to obtain 
good Tc. There is some experimental evidence for a tendency to cluster 

The effect of clustering on Tc in Gao.92Mno,o8As was studied by a simple model. To 
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FIG. 2: Dependence (K) on . in Ga,_.M„,A= and Ga,_.Cr.,A.. Sol.d lines: computed 

from the MF T" 2€\. Dotted line: extracted from spin-dynamics simulations of Eq. ^ 
Diamonds: Tc computed from the Heisenberg Cluster Variation Method. 

characterize the configurational disorder we adopt the standard Ising formahsm, and assign 
o" = ±1 to each cation site (+1 for Mn and —1 for Ga). The random (SQS) configuration 
was constructed by searching for configurations which best approximate the ideal random 
configuration for pair correlation functions Pr^. = {(JiCXj) (and some higher-order correlation 
functions) up to some fixed distance. For a random configuration, Pr = (2x — 1)^ inde- 
pendent of R. To parameterize the clustering in a simple manner, we adopted the NN pair 
correlation function Pi as a measure of clustering. Starting from an initial SQS configura- 
tion, a simulated annealing cycle was performed by generating a set of site configurations 
with increasing Pi, corresponding to longer annealing times (For simplicity, Pn{n > 1) was 
optimized to be (2x — 1)^ for each configuration.) Jij and Tc were computed by the CVM 
and MFaI^I as a function of Pi; see Fig. |3] is rather strongly reduced with increasing 
Pi. This is perhaps not surprising since increased clustering implies more distant average 
separation between atoms, which is deleterious to links in the percolation path. Even within 
the MFA Tc changes slightly, albeit for a different reason. In that case, there is an increase 
in NN pairs, which would increase Tc, but at the same time there is some increase in the 
likelihood of three- and higher body neighbors. The presence of a third neighbor has the 
effect of reducing the pairwise Jij [8], and is the origin of the factor-of-three variations in 
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the NN J in Fig. d 

We also considered the ordered by putting 1 Mn in a 24-atom unit cell, correspond- 
ing to x=0.083. In this case Pi decreases to 2/3, and Tc increases to 350K (see Fig. EI). 
Thus we conclude that ordering increases Tc, while clustering decreases T^. Perhaps not 
suprisingly, the MFA approaches the CVM result in the ordered case, since percolation 
is less critical. 
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FIG. 3: Dependence of Tc (K) on the pair correlation function Pi in Gao.92Mno.08As. The random 
(SQS) configuration corresponds to Pi=0.7056. (Two SQS structures were calculated.) Diamonds 



show Tc computed with CVM; circles show 
Pi = 2/3 corresponds to the ordered compound. 
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and triangles show T^^^^. The point at 



To conclude, we have shown that ferromagnetism is very sensitive to configurational 
disorder in DMS alloys, and that with proper treatment of disorder Tc is reasonably predicted 
by the LSDA for Gai_^.Mn^As, without needing to invoke compensating defects. We briefly 
consider two important sources of error from elements missing in the theory. First, spin-orbit 
coupling strongly reduces Tc in k-p models. We estimated its effect by computing the change 
in Tc^^^ when the L ■ S coupling is added to the LSDA hamiltonion. For Gao.92Mno.ogAs, 

TT^MFA 



Tc was reduced by ~10%. Finally, the LSDA itself will overestimate Tc somewhat 
In a future work we will present a reliable parameter-free theory that corrects the principal 
errors in LSDA — most importantly the Mn d character at Ep — and quantify the extent 
to which the LSDA overestimates T^. Finally, we conclude that the present calculations 
represent a rather strict upper bound to Tc, and that for random or clustered Gai-^Mn^jAs 
alloys, Tc > 250 K is unlikely. 

This work was supported by the Office of Naval Research. 
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